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ABSTRACT 

Rotation-powered radio pulsars are generally observed to pulse regularly in the radio 
band, but this is not the case for so-called rotating radio transients (RRATs) which 
emit only sporadic bursts of radio emission. We demonstrate that the high-magnetic 
field pulsar J1119-6127 exhibits three different types of behaviour in the radio band. 
Trailing the "normal" profile peak there is an "intermittent" peak and these compo- 
nents are flanked by two additional components showing very erratic "RRAT-like" 
emission. Both the intermittent and RRAT-like events are extremely rare and are 
preceded by a large amplitude glitch in the spin-down parameters. The post-glitch 
relaxation occurs on two different timescales (^^20 and ~210 days) and the post-glitch 
spin-down rate is smaller than the pre-glitch rate. This type of relaxation is also seen 
in an earlier, smaller glitch and is very unusual for the pulsar population as a whole, 
but is observed in the glitch recovery of a RRAT. The abnormal emission behaviour 
in PSR J1119-6127 was observed up to three months after the epoch of the large 
glitch, suggestive of changes in the magnetospheric conditions during the fast part of 
the recovery process. We argue that both the anomalous recoveries and the emission 
changes could be related to reconfigurations of the magnetic field. Apart from the 
glitches, the spin-down of PSR J1119-6127 is relatively stable, allowing us to refine 
the measurement of the braking index {n = 2.684 ± 0.002) using more than 12 years 
of timing data. 

The properties of this pulsar are discussed in light of the growing evidence that 
RRATs do not form a distinct class of pulsar, but rather are a combination of different 
extreme emission types seen in other neutron stars. Different sub-classes of the RRATs 
can potentially be separated by calculating the lower limit on the modulation index of 
their emission. Unlike other quantities, this parameter is independent of observation 
duration allowing a direct comparison with other emission phenomenon. We speculate 
that if the abnormal behaviour in PSR J1119-6127 is indeed glitch induced then there 
might exist a population of neutron stars which only become visible in the radio band 
for a short duration in the immediate aftermath of glitch activity. These neutron stars 
will be visible in the radio band as sources that only emit some clustered pulses every 
so many years. 

Key words: pulsars: individual: J1119-6127 — pulsars: general — radiation mecha- 
nisms — polarization 



1 INTRODUCTION 

PSR J 1119-6127 is a youn g, isolated radio pulsar discov- 
ered bv lCamilo etall (|2000l ) with the Parkes radio telescope 
in Australia. Its period (P = 0.41 s) has a typical value ob- 
served for pulsars, but its period derivative {P — 4.0x 10~^^) 
is among the highest known. Therefore the inferred surface 
(dipole) magnetic field strength is very large {Bs = 4.1 x 10^'^ 



G), making this pulsar one of the "high-magnetic field pul- 
sars". This field strength is comparable to the so-called 
"quantum critical field strength" of Ba = 4.4 x 10^^ g|3 
above which photon spli tting may prevent pair production 
(|Baring fc Hardin3ll998l ). Because this pulsar does not ro- 
tate particularly rapidly, the magnetic field strength at the 
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^ The surface magnetic field strength of PSR J1119— 6127 at its 
magnetic pole might actually be larger than Bct, because the 
dipole field there is twice as strong as the equatorial field Bg. 
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light cylinder (Blc = 5.7 x 10^ G) is high, but not as ex- 
treme as that of other pulsars. From its spin parameters one 
can further infer its characteristic age (tc = 1.6 kyr), which 
makes its association with supernova remnant G292.2-0.5 
likely (|Crawford et al.ll200ll ). 

The high spin-down energy loss rate {E = 2.3 x 10^® 
erg s~ ) makes this source a prom ising target for the Fermi 
7-ray satellite (jSmith et al.l |200^ V Hence it is one of the 
energetic pulsars which are currently monitored using the 
Parkes radio tel escope allowing the det ection of possible 7- 
ray pulsations (|Weltevrede et all 120101 ). During these (on- 
going) observations we noted that this pulsar has striking 
similarities with a group of neutron stars called r otating ra- 
dio transients (RRATs: [McLaughlin et al.llioo^ 'l and with 
the inter mittent pulsars (PS R B1931-I-24 being the first dis- 
covered; I Kramer et al.ll2006l). This pulsar can therefore be 
seen, like PSR B0656-l-14 (| Weltev rede et al. 200d), as an ob- 
ject that links different "classes" of neutron stars. We briefly 
introduce these different types of objects below. 

RRATs are characterized by the sporadic nature of 
their emission. Typically, detectable radio emission is ob- 
servable for less than one second per day, causing standard 
periodicity searches to fail in detecting a periodic signal. 
Therefore one has to rely on the greatest common divi- 
sor of the time between bursts to derive the underlying 
rotational period of the star. T here are about forty re- 



ported RRATs in the literature (iMcLaughlin et al 



Hessels et all l2008l: iDeneva et al.l l2009l : iKeane et al 



20061: 



2OIOI : 



Burke-Spolaor fc Bailesll2010l '). Manv RRATs have relatively 



long periods up to 7 s, suggesting these sources may be re- 
lated to the radio- quiet X-ray populations of n eutron stars, 
such as magnet ars (Woods fc Thompsonll2006l ) and isolated 
neutron stars (|Kaplanll2008l ). which are observed to have 
similar periods. The size of the RRAT population is thought 
to be several times larg er than the radio pulsar population 
iMcLaughlin et al.i2006l ). leading to inconsistencies with the 
observed supernova rate unless the re is an evolutionary lin k 
between the neutron star classes (|Keane fc Kramer] l200i ). 
Different ideas are put forward to explain the sporadic na- 
ture of detectable pulses, such as intermittent particle pre- 
cipitation towards the star from a radiation belt similar to 
those in planetary magnetospheres ||Luo fc Melrose|[2007t ) or 
re-activation of the usually inactive vac uum gap due to in 
falling circumstellar asteroidal material jCordes fc ShannonI 
[2008^ 

IWeltevrede et al.l (|2006l ') argue that PSR B0656+14 
would have been classified as a RRAT were it not one of 
the most nearby pulsars. One is therefore left with the ques- 
tion if PSR B0656+14 should be called a RRAT? If not, 
it implies that at least some of the "RRATs" are physi- 
cally very similar to PSR B0656+14, except that they are 
more distant. In any case, it is clear that defining RRATs to 
be pulsars that are not detectable via periodicity searches 
is insufficient and that one has to define the different neu- 
tron star classes b ased on physical properties instead (e.g. 
iKeane et aklboid '). 

Based on the similarities between PSR B0656-I-14 and 
the RRATs. IWeltevrede et~al] jiooe ) predicted that at least 
some of the RRATs would show additional emission which 
is much weaker than their bright individual pulses and 
that when the individual bursts of RRATs are averaged 
they will form a "profile" much broader than the individ- 



ual bursts. Indeed, for instance the RRAT PSR J0627-I-16, 
which was found via its individual pulses, showed weak emis- 
sion in a follow-up observation (|Deneva et al.l |2009| ) and 
the sum of the individual bursts of the RRAT PSR J1819- 
1458 produces a triple peaked p ulse profile which is muc h 
wider than the individual bursts (|Karastergiou et al.ll200ij ). 
The latter RRAT shows a X-ray spectrum which is con- 
sistent with thermal e mission from a cooling neutron star 
( Reynolds et al. 200^ : iGaensler et al. I I2OO7.) and, like PSR 
B065 6-I-14, it shows X-ray pulsations ( McLaughlin et al] 
|2003). As will be discussed later in this paper, the pecu- 
liar emission properties of PSRs J1119-6127 and B0656-f 14 
might be linked to the erratic emission seen for other young 
and energ etic pulsars, such as the Vela pulsar and PSR 
B 1706-44 (| Johnston et al.ll200ll : Dohnston fc Romanill2002l '). 
thereby providing some generalization between the different 
phenomenon. 

The final class of neutron stars we will discuss in the 
introduction are the intermittent pu lsars. The archetype is 
PSR B1931-I-24 (|Kramer et al.ll2006l ). which is only active 
for a few days between periods of roughly a month during 
which the pulsar is not detectable. These sudden switches 
are accompanied by changes in the spin-down rate of the 
neutron star rotation, suggestive of significant changes in 
the torque generated by magnetospheric currents. It now 
has become clear that many more objects show pulse pro- 
file changes concu rrent with changes in the spin-down rate 
(|Lvne et all bOld l and it has been suggested that all the 
so-called timing- noise could be ascribed to this effect. In- 
deed the spin parameters of PSR B1931-I-24 (P = 0.81 s 
and P = 8.1 x 10~^*) and other physical quantities such as 
Bs = 2.6 x 10^^ G,E = 5.9 x 10^^ ergs"^ and = 1.6 Myr 
are not remarkable, supporting the idea that the correlated 
profile and spin-down switches could be very common in the 
pulsar population. It is therefore natural to extend this idea 
to all pulsars that switch off (or change their profile) on a 
regular basis during so-called "nulls" or "mode changes". 
The nulls of some pulsars are very short (a period or less), 
while oth ers are only switc hed on for a few percent of the 
time fe.g. lWang et al.ll2007l '). 

In the next sections we give some details about the ob- 
servations (Sect. ^ and we will describe a transient profile 
change in Sect. [31 This is followed by the analysis of the 
individual pulses (Sect. [4]), polarization (Sect. [5]) and pulsar 
timing behaviour (Sect.[S]). The paper ends with a discus- 
sion of the properties of PSR J1119-6127 in the light of the 
different neutron star classes (Sect. [7| and the conclusions 
(Sect.iJ. 



2 OBSERVATIONS 

All the discussed observations were carried-out using the 64- 
m Parkes radio telescope in Australia and the majority of 
the data wer e obtained as part of th e aforementioned tim- 
mg program (IWeltevrede et al.ll2010l ). This program started 
in April 2007. Each pulsar is typically observed once per 
month at 20 cm and twice per year at 10 and 50 cm, al- 
though this pulsar is only detected at 10 and 20 cm. The 
used receivers were the centre beam of the 20 cm multi- 
beam receiver (which has a bandwidth of 256 MHz and has 
a noise equivalent fiux density of ~35 Jy on a cold sky) 
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Figure 1. The only observation known to us showing PSR J1119— 6127 in its double-peaked mode. In both panels the intensity is shown 
in grayscale. Left-hand panel: Both profile components are identically affected by dispersion (an observing frequency dependent delay 
caused by the interstellar medium), resulting in the two parallel diagonal bands in this panel. Right-hand panel: The double-peaked pulse 
profile was visible during the entire 3 minute observation. 



and the 10/50 cm receiver (which has at 10 cm a band- 
width of 1024 MHz and has a noise equivalent flux density 
of ~49 Jy on a cold sky). The data were on-line folded at 
the pulse period by a digital filterbank and were dumped 
every 30 seconds on hard disk. For details about the polar- 
ization calibration procedure and the used method to sum all 
the individual observati ons into a high signal-to-noise " stan- 
dard" profile we refer to lWeltevrede fc JohnstonI (|200 8l). For 
the majority of the observations the analogue filterbank was 
used to record the pulsar signal in parallel with a time res- 
olution of 250 ^s. Although only total intensity is recorded 
it has the advantage that the data are not folded on-line, 
hence allowing the analysis of individual pulses. 

In addition we made use of archival Parkes data 
recorded between 1998 February 7 and 2007 February 12. 
For de tails about these observations we refer to lCamilo et al.l 
l|2000h in which some of these analogue filterbank data are 
published. Together with data from our timing program up 
to 2010 June 16, this results in a timing history of more than 
12 years. 



3 THE TRANSIENT PROFILE COMPONENT 

The pulse profile of PSR J1119- 6127 is weU known 
to be single-peake d (se e e.g. ICamilo et al.l I2OO0I : 
Ijohnston fc Weisberd I2OO6I ). The observation recorded 
in 2007 June 17 at 02:35:19 Universal Time therefore 
immediately triggered our interest in this object as it 
shows a clear double-peaked profile (see Fig. [T]). Despite 
years of trying we never observed a similar event again. 
Hence, it was important to convince ourselves that the 
signal was pulsar related rather than caused by something 
instrumental or by terrestrial radio frequency interference 
(RFI). As one can see in the left panel of Fig.[T] both peaks 
of the pulse profile show an identical frequency dependent 
delay (the gray horizontal bands are frequency channels 
affected by RFI and were excluded from the analysis). This 



frequency dependent delay is exactly the same as what 
is normally observed for this pulsar and can be perfectly 
explained in terms of dispersion by the interstellar medium. 
This, combined with the fact that the double-peaked 
structure is seen during the whole observation after folding 
at the known pulse period (see right panel of Fig. [1} rules 
out RFI as the cause of the additional component. The 
signal was recorded using two independent digital filterbank 
backendfl and the recordings look identical as expected. 
All observations before and after this particular observation 
(of different pulsars) used an identical system set-up and 
do not show any sign of a similar effect. Therefore the 
appearance of the additional component is clearly pulsar 
related. 

With only one instance of the pulsar showing its double- 
peaked mode and no observation of an actual transition it 
is impossible to reliably estimate the occurrence rate of this 
transient phenomenon. The pulsar was double peaked during 
the whole 3 minutes of the observation and the observations 
49 days earlier and 35 days later look normal, implying a to- 
tal duration of the event of at least minutes, but shorter than 
months. In order to maximize the chance of catching the pul- 
sar again during a similar event we substantially increased 
the observing duration (up to half an hour) in subsequent 
observing sessions. In total we have observed this pulsar 76 
times in between 2007 April 30 and 2010 June 16, typically 
twice a day every month with a total duration of 7.9 hours. 
In addition we inspected archival Parkes data recorded be- 
tween 1998 February 7 and 2007 February 12. None of these 
256 additional observations (with a total duration of 41.9 
hours) shows a similar event. Therefore the double-peaked 
mode is clearly very rare as it is only seen in 0.10% of the 
inspected data. 



The used backends are known as "DFBl" and "DFB2". 
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Figure 2. Top panel: Averaged pulse profile of a single 22 min- 
utes long observation at a wavelength of 20 cm. The time res- 
olution is reduced to 200 bins across the pulse period and the 
peak of the profile is put at pulse longitude 0°. Bottom panel: 
The strong (non-consecutive) individual pulses observed during 
the same observation. Notice that at the pulse longitude of the 
bright individual pulses there is no integrated power visible, which 
is reminiscent of the emission of RRATs. 



4 SINGLE-PULSE ANALYSIS: THE 
RRAT-LIKE COMPONENTS 

4.1 20-cm observations 

At first glance this pulsar is not a promising source to anal- 
yse its individual pulses because it is quite weak (the aver- 
age single-pulse signal-to- noise is less than 1). Nevertheless, 
careful investigation by eye of all the observations for which 
we had single-pulse recordings available revealed a handful 
of strong pulses. In total only four pulses at an observing 
wavelength of 20 cm were strong enough to be clearly de- 
tectable (see Fig. [2]). They all occurred during a 22 minutes 
long observation recorded at 2007 August 20, starting at 
22:11:05 UT. 

Strong pulses were visible only in one of the 217 observa- 
tions (with a total duration of 28.3 hours) for which we have 
single-pulse data (and which were not too badly affected by 
RFI). In other words, only during 4 of the 246,892 inspected 
neutron star rotations the intensity of the radio beam was 
strong enough to be detectable with the Parkes telescope. 




5 10 15 

Frequency [Hz] 



20 



a ° 



5 10 15 

Frequency [Hz] 



20 



Figure 3. Top panel: Spectrum of the time-series of the obser- 
vation containing the four strong pulses shown in Fig. [2] The 
periodicity of the pulsar signal clearly stands out as a spike at 
2.45 Hz and spikes at higher frequencies harmonically related to 
the fundamental frequency. Bottom panel: The spectrum of the 
same time-series after replacing the signal at the pulse longitudes 
corresponding to the normal emission by Gaussian noise thereby 
only keeping the strong pulses. None of the visible spikes corre- 
spond to the periodicity of the pulsar. The normalization of the 
spectral power is the same in both panels. The decrease in spec- 
tral power below ~ 1 Hz is due to a software-filter which removes 
baseline variations. 



Assuming a constant event rate for the occurrence of strong 
pulses, Poisson statistics predicts an extremely low proba- 
bility of 3 X 10"^ that four (or more) pulses occur during 
a single 3211 pulse period long observation by pure chance. 
This strongly suggests that the event rate must have been 
much higher during that particular observation. 

Looking at Fig. [2] one can see that the strong individ- 
ual pulses all occur ~ 35° in pulse longitude earlier than 
the peak of the normal profile. Surprisingly, at the pulse 
longitudes one can detect individual pulses, the integrated 
profile does not show significant emission. This is a defining 
characteristic of RRATs, as it implies that this kind of emis- 
sion cannot be detected using standard periodicity searches. 
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To prove this point, we took the observation with the four 
strong pulses and replaced the emission within the pulse 
longitude range corresponding to the normal emission with 
Gaussian noise with a root mean square (RMS) equal to the 
off-pulse RMS. The spectrum based on this time-series is 
shown in the bottom panel of Fig. [3l while the spectrum in- 
cluding the normal emission is shown in the top panel. One 
can clearly see that the fundamental frequency (2.45 Hz) 
and the higher harmonics are only visible when the normal 
emission is included. Despite the strongest pulses are ob- 
served in the leading component, the underlying periodicity 
(or the underlying emission itself) is not detectable at those 
pulse longitudes. Therefore the leading component can be 
qualified to have RRAT-like emission, although as will be 
discussed later, there are also important differences. 

The strongest of the four detected pulses is 23 times 
brighter than the average pulse of the normal emission dur- 
ing the same observation. Assuming the normal emission 
of the pulsar was not abnormally bright or weak during 
that particular observation (and there is no reason to believe 
so based on the signal-to-noise ratio of the pulse profile of 
that observation compared to others), this corresponds to 
18 mjy averaged out over a whole pulse period (using the 
published mean flux d ensity of 0.80 mJy at 1400 MHz by 
[Manchester et al]|200ll ) or a peak flux density of 0.6 Jy us- 
ing the full single-pulse width of 11°. 

To investigate how "normal" the radio emission in 
the main-peak following the RRAT-like component of PSR 
J1119-6127 is, we attempted to calculate the modulation 
index for this pulsar. The modulation index is independent 
of the observing system and quantifies the broadness of an 
amplitude distribution. It is defined to be 
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where ai is the standard deviation of the observed pulse 
intensities, (/) is the average pulse intensity and A^tot is the 
total number of stellar rotations observed. This quantity can 
be calculated as function of pulse longitude by considering 
the intensities at a specific longitude rather than the pulse 
intensities integrated over pulse longitude. We have chosen 
to calculate the modulation i ndex in the Fourier domain 
instead, which as explained bv lEdwar ds fc St appersI (|2002l ') 
allows a correction for interstellar scintillation. 

It turned out that the signal-to-noise ratio in an indi- 
vidual observation is not high enough to detect a signifi- 
cant modulation index. To circumvent this problem we have 
added together all the individual pulse recording at a wave- 
length of 20 cm. In order to do this the average off-pulse 
intensity of each observation was subtracted, the data were 
re-sampled to have an equal number of bins across the pulse 
period and the intensities were scaled to make the off-pulse 
RMS equal to 1 (thereby assuming that the system temper- 
ature was the same for all observations). The individual ob- 
servations were then aligned by cross-correlating their pulse 
profiles with an anal ytic template based on von Mises func- 
tions (see e.g. Welte vrede fc Johnston 2008i). The resulting 
pulse stack is based on more than 28 hours of data and shows 
a modulation index of m ~ 0.5 in the centre of the profile, 
flaring up at the edges. It is likely that the modulation index 
is slightly lower in reality because the addition process can 
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Figure 4. Top panel: Pulse profile of a single 6 minute obser- 
vation at a wavelength of 10 cm. The time resolution is reduced 
to 300 bins across tlie pulse period. The peak of the pulse pro- 
file is placed at zero longitude. Bottom panel: The strong (non- 
consecutive) individual pulses observed in the same six minute 
observation are plotted in grayscale. 



add some artiflcial variations to the signal. Nevertheless, the 
modulation index of PSR J1119-6127 is e ntirely consistent 
to th ose observed for most pulsars (e.g. i Weltevrede et al.l 
l2006h . 



4.2 10-cm observations 

PSR Jl 119-6127 has only been observed sporadically at a 
wavelength of 10 cm during our timing program and at 
a wavelength of 50 cm the pulsar is too weak to be de- 
tectable for the used integration times. It must be noted 
that the combination of a high dispersion measure, worse 
system temperature, smaller band width and a much higher 
sky temperature might well be preventing us from detect- 
ing the pulsar at 50 cm. Because the individual pulses were 
not always recorded at those wavelengths, there are only 
two individual pulse recordings available at 10 cm. One of 
these recordings, the 6 minute long 2007 July 23 observation. 
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Figure 5. The high signal-to-noise standard profile showing the pulsar in its single-peaked mode (left panel) and the profile in its 
transient double-peaked mode (middle panel). The top plots show total intensity (solid line), linear polarization (dashed line) and 
circular polarization (dotted line). Zero longitude corresponds to the measured location of the steepest gradient of the PA curve (shown 
in the bottom panels). The bottom panel of the middle panel shows the RVM fit through the PA points. The exact same curve is shown 
in the left-hand panel, without any additional fitting applied. The perfect match, including the absence of any offset between the data 
and the model, unambiguously shows that it is the leading component of the double profile which corresponds to the single-peaked mode. 
Right panel: The total intensity profiles of the single-peaked mode (solid line) and the double-peaked mode (dotted line), as well as their 
PA points, are overlaid. The bottom panel illustrates again the perfect match of the PA-points. 



showed a number of clearly detected bright individual pulses 
(bottom panel Fig. 2)), while a 4 minute long 2009 August 9 
observation does not show any sign of similar pulses. 

The strong pulses appear to occur at three distinct 
pulse longitudes (approximately at —30°, —10° and 50°). 
Notice that the pulse profile of this particular observation 
shows structure at these pulse longitudes associated with the 
strong individual pulses (top panel of Fig. [4} . This suggests 
that the strong pulses appear before and after the main peak 
of the observation. The relative alignment of this profile with 
respect to the single-peaked and double-peaked profile ob- 
served at a wavelength of 20 cm will be discussed in the next 
section in more detail. In total there are eight observations 
at a wavelength of 10 cm for which we have an integrated 
pulse profile and this observation is the only one which shows 
evidence for additional components in the pulse profile. The 
fact that the strong pulses do have corresponding peaks in 
the pulse profile suggests that at this frequency (or possibly 
epoch) the strong pulses are less RRAT-like than the strong 
pulses observed at a wavelength of 20 cm. 

Independent of whether single pulses occurred in one 
of the two, or eight, of the 10-cm observations, given how 
rare the detection of individual pulses is at a wavelength of 
20 cm and the fact that at 10 cm more individual pulses 
are detected suggests that strong pulses occur more often 
at higher frequencies. This is consistent with the fact that 
there are no individual pulses detected in two 3 minute ob- 
servations made at a wavelength of 20 cm one day earlier 
and two days later. Only if the event rate is high for periods 
of time less than a day and we have been lucky during the 



observation at a wavelength of 10 cm the event rates could 
be the same at the two frequencies. 



5 POLARIZATION PROFILES 

It is not immediately obvious which of the two compo- 
nents of the double-peaked profile corresponds to the single- 
peaked profile and therefore it is not clear if and how the 
newly emerged component relates to the RRAT-like emis- 
sion. In principle one could use pulsar timing methods to 
resolve this issue, but, as will be discussed in Sect. [6] a large 
amplitude glitch makes this impossible. We therefore have 
to resort to other measurements to find out how to time 
align the normal single-peaked pulse profile and the double- 
peaked profile. We will show in the this section that the 
position angle (PA) of the linear polarization can be used 
to find this relative alignment. Polarization data are further 
used to constrain the viewing geometry and the radio emis- 
sion height. 

5.1 Alignment of the profile components 

Virtually all pulsars with E > 2 x 10^^ ergs~^ have a lin- 
ear p olarization fraction over 50% jWeltevrede &i JohnstonI 
l2008l : see also e.g. lvon Hoensbroech at alj|l998l ). In addition, 
most of these pulsars show a PA which varies smoothly as 
a function of pulse longitude forming a S-shape. This pul- 
sar is no exception as one can see in the left panel of Fig. 
[5] The middle panel shows the double-peaked profile. No- 
tice that the PA values of the pulsar in its normal mode 
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clearly match those of the leading peak of the double-peaked 
profile, suggesting that the trailing peak is the component 
which has appeared in the double-peaked mode. Indeed a 
cross-correlation between the two profiles shows the highest 
degree of correlation for such an alignment. This is caused 
by the slight asymmetry seen for the profiles. In the single- 
peaked mode the trailing side of the profile is slightly steeper 
caused by a weak emission component overlapping with the 
leading edge of the main peak (visible at —35° pulse longi- 
tude). The leading component in the double-peaked mode 
has also a shallower leading edge, although in that case there 
is no evidence for a weak component being responsible. 

From the fact that the PA curve is very similar in the 
single- and double-peaked mode one can conclude that the 
PA provides a much stronger foundation to base the align- 
ment of the two profiles on than the total intensity. In or- 
der to quantify the relative alignment of the two profiles we 
have fitted the data using the rotating vector model (RVM; 
iRadhakrishnan fc Cookelll969l ). In this model the PA of the 
linear polarization ij) depends on the pulse longitude (j) via 



tan(V' - V'o) = 



sm a sm 



sin C, cos a — cos C, sin a cos ( 



(2) 



where Vo and (jio are the PA and pulse longitude at the 
location where the ^(0) has its steepest gradient (and an 
infiection point). In addition, the shape of the PA curve 
depends on the emission geometry of the pulsar, which is 
described by two angles. Firstly the angle a between the 
magnetic axis and the rotation axis and secondly the angle 
C, between the line-of-sight and the rotation axis. A related 
angle is the impact parameter P = C, — a, which is the angle 
between the line-of-sight and the magnetic axis at its closest 
approach. 

According to the RVM, corresponds to the pulse lon- 
gitude at which the line-of-sight intersects the fiducial plane, 
the plane containing the rotation and magnetic axis. There is 
a longitude-shift one has to take into account if the emission 
height hem is not negligible compared to the light cylinder 
distance, as will be discussed later. The longitude (^o could 
act as a natural zero point for the pulse longitude scale if 
one uses the PA to align different profiles. This is achieved 
by determining <^o by fitting Eq. [2] through the PA points. 
The profiles are then rotated such that the <j)o corresponds 
to 0° in the plots of Fig. [S] Note that the alignment is in- 
dependent of the actual value of the PA, it only depends on 
the pulse longitude corresponding with the steepest gradi- 
ent of the RVM curve. This allows a self-consistency check of 
the determined alignment by over-plotting the RVM curve 
obtained from the fit of the double-peaked profile on the 
PA points of the single-peaked profile. One can see in Fig. 
[5] that this RVM curve matches both PA curves perfectly, 
something which is also evident from the right-hand panel 
which shows the PA points of the two profiles overplotted. 
This shows that the measured position of the steepest gradi- 
ent of the PA curve is robust, because otherwise there would 
be a vertical (or horizontal) offset between the PA points of 
the two profiles. 

Now that the single-peaked profile is aligned with the 
double-peaked profile, the leading RRAT-like component 
can be placed on the same longitude scale by aligning the 
main peaks. This is shown in Fig. [6] which shows both the 
sum of the four strong pulses observed at 20 cm (dotted line) 



c 
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Figure 6. Plot showing the profiles obtained by averaging the 
four strong pulses shown in Fig.[2](dotted line), the double-peaked 
profile (thick gray line) and the abnormal profile at 10 cm (solid 
line, with a vertical ofl'set from the other two profiles). The pulse 
longitude scale is identical to that in Fig. [5] and the method for 
aligning the 10-cm profile is explained in the text. The peak in- 
tensities of the profiles are normalized. 



and the double-peaked profile (thick gray line) at the same 
wavelength. The strong pulses represent a distinct profile 
component with a pulse longitude separation from the main 
peak roughly equal to that between the main peak and the 
transient trailing peak of the double-peaked profile. 

The remaining question is how the components seen in 
the 10-cm observation discussed in Sect. |4]2] fit in this pic- 
ture. The profile is far too weak to find the pulse longitude 
at which the gradient of the PA curve is steepest. Never- 
theless, the value of the PA can be used to roughly align 
the profile to the PA curve observed at 20 cm. It follows 
that the main peak of the 10-cm profile overlaps roughly 
with the single-peaked profile observed at 20 cm. In Fig. [6] 
the 10-cm profile is shown together with the double-peaked 
profile and the RRAT-like component seen at 20 cm. This 
alignment suggests that the bright pulses at the far leading 
edge of the 10-cm profile correspond with the RRAT-like 
component seen at 20 cm (dotted line). The strong pulses 
seen at 10 cm at the far trailing edge of the profile appear to 
form another profile component, again with a similar compo- 
nent separation as between the other components. It should 
be emphasized that because of the low signal-to-noise ratio 
of the 10-cm profile, its alignment with the other profiles 
should only be considered as a rough indication. 

We therefore believe that although the pulse profile of 
PSR J1119-6127 is most of the time single-peaked, there 
are in fact at least four components. In addition to the nor- 
mal emission, the profile has shown once a transient trailing 
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Table 1. This table summarizes the considered choices for the 
location of the magnetic axis with respect to the pulse profiles of 
PSR J1119— 6127. The second and third column denotes those lon- 
gitudes in the single- and double-peaked mode respectively using 
the longitude shown used in Fig. [5] The next three columns indi- 
cate the implied emission height difference between the two modes 
and the emission heights themselves. The last column states which 
component of the double-peaked profile is the new component. In 
case II the single-peaked profile splits into two components. Both 
cases I and IV are argued to be plausible, case IV being our pre- 
ferred solution. 



Case 


0sgl 


</>dbl 


A/idbi 


^sgl 


^dbl 


New comp. 








[km] 


[km] 


[km] 




I 


-16° 


-21° 


400 


1400^ 


1800 


Trailing 


II 


-16° 


-6° 


-900 


1400"~ 


500 


Split 


III 


-16° 


+7° 


-2000 


1400 


-600 


Leading 


IV 


-6° 


-6° 





500 


500 


Trailing 


" An 


emission 


heif 


;ht of 


2200 km 


was 


derived by 



with the steepest 



radient lag- 
ging the peak of the total intensity profile by 26°. We attribute 
the discrepancy between their result and ours to the lower quality 
of their data and note that their pulse profile and PA curve are 
fully consistent with those presented in this paper. An emission 
height of 2100 km was derived bv I Weltevrede fc JohnstonI l l2008i) , 
however in their table all emission heights derived from PA curves 
are overestimated by a factor 2 because of a miscalculation. 

peak. These two peaks are flanked by two RRAT-like com- 
ponents which are usually inactive. 



5.2 Emission heights 

5.2.1 20-cm emission 

The single-peaked proflle does not perfectly match the lead- 
ing component of the double-peaked profile as shown in the 
right-hand panel of Fig. [S] There is no a priori reason why 
they should match better, given that the shapes of the pro- 
files are slightly different to start off with. On the other hand, 
there are reasons for the alignment of the profiles not to be 
perfect. If the emission height is different in the two states 
of the pulsar, which potentially could be the root-cause of 
the different profile shapes, one does expect a deviation in 
the alignment of the profiles. This is because only for a zero 
emission height 4>o corresponds to the pulse longitude at 
which the line of sight is passing the fiducial plane. For finite 
emission heights co-rotation of the emitting region makes the 
infiection point of the PA curve ((^o) to be delayed with re- 
spect to the pulse profile. The relative shift A</) between the 
pulse profile and the inflection point of the PA curve re lates 
to the emission height /lem via (|Blaskiewicz et ai]|l99ll ) 



1 Pc 

ftom = -jRhcA<j} = —A(j). 

4 Stt 



(3) 



Here P is the spin period of the pulsar, c is the speed of 
light and _Rlc is the Hght cylinder radius (19,500 km for 
this pulsar). Therefore, if the emission heights are different 
in the single- and double-peaked modes, the pulse proflles 
will not be perfectly aligned when using the inflection point 
of the PA as the reference point. 

In order for the centroids to overlap, the normal pro- 
flle should be shifted by 5° to earlier longitudes compared 



to the double-peaked proflle. This corresponds to an emis- 
sion height which is 400 km (2% of the Hght cylinder ra- 
dius) lower in the single-peaked mode (case I in Table [T]). 
Such an emission height difference would make both the PA- 
points and the total intensity proflle overlap. If one allows 
the emission height difference between the two modes to 
be bigger, more extreme offsets could be considered. For in- 
stance, one could envision the flducial plane to correspond to 
the centre of both the single-peaked and the double-peaked 
proflle, thereby making both proflles symmetric around the 
magnetic axis (case II). Or even more extreme, maybe the 
single peak corresponds to the trailing component in the 
double-peaked mode (case III). In the flnal case we con- 
sider the magnetic axis corresponds to the symmetry point 
of the double-peaked proflle without a differential emission 
height (case IV). With the made assumptions about which 
pulse longitude of the pulse proflle corresponds to the fldu- 
cial plane, Eq. [3]can also be used to constrain the emission 
height itself rather than just height differences. The result- 
ing emission heights of the single- and double-peaked mode 
can be found in Table [T] 

Given the fact that the component width is not very 
different in the normal mode and the double-peaked mode, 
one could argue that the emission height difference between 
the two modes must be relatively small compared to the to- 
tal emission height. For case II (i.e. the single-peaked proflle 
splits in the double-peaked mode, arguably the most sym- 
metric scenario) one would expect the overall pulse width 
to be roughly 1.7 times larger in the single-peaked mode 
because of the curvature of the dipolar magnetic fleld lines. 
This effect must be very precisely counter-balanced by hav- 
ing a smaller fraction of the polar cap being active during 
the single-peaked mode in order to have similar component 
widths, at which point the model becomes contrived. This 
problem becomes even worse in case III, which is already 
unlikely given the negative derived emission height in the 
double-peaked mode. 

In case IV a differential emission height cannot explain 
the apparent small offset between the single-peaked proflle 
and the leading component of the double-peaked proflle, as 
an essentially zero emission height would be required for 
the single-peaked mode. Nevertheless, the striking mirror 
symmetry of the double-peaked proflle can be seen as strong 
evidence for the location of the flducial plane at that pulse 
longitude. 

Therefore we conclude that, even when differential emis- 
sion heights are considered, the most likely scenario is that it 
is the trailing peak which has emerged in the double-peaked 
mode (either case I or IV in Table [T|). Possibly (although 
not necessarily) the flducial plane corresponds to the exact 
centroid of both the proflle of the single-peaked mode and 
the leading component in the double-peaked mode (case I). 
Otherwise the emission height might be identical in the two 
modes with the flducial plane corresponding to the centre of 
the double-peaked proflle (case IV). 

5.2.2 10-cm emission 

Fig. [7| shows the standard proflles obtained from the timing 
program at a wavelength of 20 and 10 cm. The longitude 
scale used for the 20-cm proflle is the same as in Fig. [5] and 
the alignment of the 10-cm proflle is obtained by matching 
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Figure 8. The map resulting from fitting the RVM model to 
the polarization data is shown in grayscale. The solid contours 
indicate where the 'x^ is two, three and four times larger than 
the minimum value. Overlaid are the p contours assuming the 
double-peaked pulse profile fills the last open field line region. 



Figure 7. The normal profile of PSR Jllf 9-6f 27 at a wavelength 
of 20 cm (solid line) and at 10 cm (dotted line). The steepest gra- 
dient of the PA-swing measured at a wavelength of 20 cm occurs 
at 0° longitude. The PA-points of the 10 cm data were aligned by 
eye to those at 20 cm and the PA-points at both frequencies are 
overlaid in the bottom panel. The alignment of the profiles at the 
two frequencies is, within the uncertainties given the relatively 
low signal-to-noise ratios, the same as in Fig. [O] 

the PA curves of the 10- and 20-cm profiles. The centroid 
of the 10-cm profile appears to be shifted to later longi- 
tudes compared to the 20-cm profile by about 3° . Using the 
same methodology as described earlier, this suggests that 
the emission height is 300 km lower at a wavelength of 10 
cm. This height diflFerence implies a difference in the ex- 
pected half opening angle p of the beam at the two wave- 
lengths when using the equation 

In this equation p is related to the opening angle of the 
last open dipole field lines at the emission height (e.g. 
iLorimer &: Krameiil2005l ) , thereby assuming that the full po- 
lar cap is active. 

The difference in emission height implies that the emis- 
sion height at a wavelength of 20 cm is a factor of 1.2 and 
2.0 higher than that of the 10-cm emission for cases I and 
IV of Table [T] respectively, which translates in a factor 1.1 
and 1.4 difference in p. The latter factor is similar to the 
observed difference in pulse width at the two wavelengths, 
therefore supporting scenario IV unless the line of sight is 
grazing the edge of the emission cone. It should be stressed 
that this is not a very direct argument, as it involves a chain 
of assumptions such as the assumption that the same field 
lines are producing radio emission at the two frequencies. 
Nevertheless, the data are at least consistent with case IV. 

5.3 Viewing geometry 

The orientation of the radio beam with respect to the line 
of sight and the rotation axis can be quantified with the 



angles a and C,. These angles can be constrained first of 
all by fitting the RVM curve (Eq. ^ . This is visualized in 
Fig. [5] which shows the of the RVM fit. As is often the 
case, the contours are banana-shaped, caused by a lack in 
pulse longitude coverage of the PA points. Therefore, not 
surprisingly, the best constraint is obtained by using the PA 
curve of the double-peaked profile rather than the single- 
peaked standard profile despite the much higher signal-to- 
noise ratio of the latter. 

An additional constraint follows from estimated half 
opening angle p of the beam (Eq. |4} , which constrains the 
angles a and C, because p is related to the pulse width W 
via 

cos p — cos a cos ^ + sin a sin cos {W/2) , (5) 

assu ming the radio beam is symmetric about the magnetic 
axis (|Gil et al.lll98i : iLorimer fc Kramedl2005l ). The dotted 
contours in Fig. |8] indicate the p values as derived using 
Eq. O using the observed pulse width of the double-peaked 
profile (which we take to be the full longitude range over 
which we observe significant emission). 

In case I the radio beam is asymmetric in the double- 
peaked mode (but symmetric in the single-peaked mode). 
Taking the outer edge of the trailing component of the 
double-peaked profile to correspond to the last open field 
line and using p = 26° (the expected opening angle for an 
emission height of 1800 km in the double-peaked mode) it 
is derived that a ~ 30° — 40°. Note that this implies that 
the strong pulses at the trailing edge of the profile observed 
at 10 cm are produced outside the polar cap region. If all 
emission occurs from within the open field line region the 
effective pulse width is larger making a smaller. 

For case IV in Table [T] the solution we preferred in 
subsection [521 P = 14° (Eq.lHl. This suggests a = 20° - 30° 
(see Fig. [8]) assuming the double-peaked profile fills the open 
field line region. If the RRAT-like emission comes from open 
field lines a more aligned geometry (a ~ 17°) is required to 
explain the larger overall pulse width (see Fig. [6} . 
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6 TIMING 

In order to characterise the spin-behaviour of a pulsar one 
needs to measure the time-of-arrival (TOA) of the radio 
pulses as a function of time. This is done by cross corre- 
lating the pulse profiles of the individual observations with 
a high signal-to-noise radio standard profile (shown in the 
left panel of Fig. [S}. Clock corrections were applied to the 
TOAs which were then converted to arrival ti mes at the solar 
system ba rycentre using the DE405 model (jEdwards et al.l 
I2OO6I: [Standish 199a ') using the TEMP02 timing package 



I Hobbs et al 



200d '). 



Next, the barycentric TOAs are compared with a timing 
model and its parameters are refined by minimising the so- 
called timing residuals. The refining process has been done 
using custom software which in essence works identical to 
timing packages such as TEMP02. The minimisation proce- 
dure makes use of the downhill simplex method and the un- 
certainties are derived from the correlation mat rix which is 
estim ated by calculating the Hessian matrix (e.g. lPress et al.l 
I1992I ). 

The most basic timing model for the rotation of a pulsar 
is a truncated Taylor series: 

(P{t) = (Po + v{t - to) + ^i^{t - tof + ii>(t - ^o)^ (6) 

where (j){t) is the rotational phase as function of time. The 
parameters (po, v, v and v are the rotational phase, spin 
frequency and its first two time derivatives, all defined at 
epoch io, usually taken to be during the time covered by the 
observations. 

PSR J1119— 6127 presents a typical timing behaviour 
of a pulsar of its age in the sense that like the Crab and 
other very young pulsars it exhibits a relatively stable spin- 
down evolution despite the significant timing noise 
observed. This allows a good estimate of its braking index as 
will be discussed later. During the time of the observations 
the rotation of PSR J1119— 6127 has been interrupted by 
two major glitches exhibiting a rather particular behaviour. 

A pulsar glitch is a sudden, normally un- 



resolved, 



increase of the rotation frequency 



IIRadhakrishnan fc Manchesteilll969l 'l. Sometimes, glitches 
are observed to be followed by an exponential-like recovery 
towards the former rotational state; phenomenon that 
has been interpreted as a signature of the presence of a 
superfluid in the interior of the star l|Bavm et al.lll96^ ). In 
many cases this recovery is either too small or the sampling 
too sparse to quantify its effects. 

The two large glitches that we report here occurred in 
2004 and 2007; the first one between MJD 53279 and 53306 
(October 1 and 28), and the second one in 2007 between 
MJD 54220 and 54268 (April 30 and June 17). All TOAs 
between MJD 50850 and 55364 were modelled using a single 
model that includes these glitches. The phase difference after 
a glitch is modelled using an additional phase </>g(t), which 
describes the deviation from Eq.[S] caused by the glitch. The 
phase after the event is parametrised by permanent jumps 
in frequency (At-p) and frequency derivative (Az>p) plus an 
exponential decay with time constant r^: 




5.2x10* 



5.3x10* 5.4x10' 
MJD [days] 



-AvdTde 



-(t-tg)/Td 



(7) 



Figure 9. Top panel: The timing residuals obtained by only mod- 
elling the permanent jumps in rotational frequency and its time 
derivative after the two glitches (their epoch is indicated by the 
dashed lines). The ordinate is expressed in turns. Middle panel: 
Timing residuals obtained by including a single exponential glitch 
recovery after each glitch. Bottom panel: Timing residuals ob- 
tained by including a glitch recovery on a second faster timescale 
after the second glitch. 



where tg corresponds to the glitch epoch and Aud is the step 
in frequency which is exponentially recovered. The step in 
rotational phase A(j} is not a physical jump in phase, but it 
corrects for a tg which is generally not specified accurate to a 
fraction of a rotational period. The instantaneous frequency 
jump at the glitch epoch is given by Avg = Avp + Aud 
and the corresponding frequency derivative jump is Avg — 
Aup + Aud, where Ai/d = -Avd/rd- 

During the minimisation process the epoch of the two 
glitches were kept fixed at the values MJD 53290 and 54240, 
respectively. The top panel of Fig. [2] shows the timing resid- 
uals after applying the above described model having set 
Avd ~ 0, i.e. neglecting any possible glitch recovery. It is 
clear that the amplitude of the residuals is much larger after 
both glitches, indicating that the glitches are poorly mod- 
elled. This, in combination with the observed pattern of the 
residuals, is suggestive of a post-glitch transient stage. Ac- 
cordingly, the residuals are considerably improved when ex- 
ponential recoveries are allowed for the two glitches (middle 
panel of Fig. Nevertheless, the residuals immediately af- 
ter the 2007 glitch show significant deviations, indicating 
that the recovery still has not been properly modelled. If a 
second exponential recovery term with a shorter time con- 
stant is included for this glitch, the residuals become much 
smaller (bottom panel of Fig. [9}. The parameters of this 
last fit can be found in Table [5] Because the recovery pro- 
cess may not yet be complete at the time of writing, the 
parameters for the 2007 glitch may change slightly if new 
data are added in the near future. 

As mentioned earlier, the timing analysis can in prin- 
ciple be used to identify which component of the double- 
peaked profile corresponds to the normal profile. This can 
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Figure 10. Top-left panel: The measured rotational frequency of PSR J1119-6127 as a function of time. The frequency is decreasing, 
as is expected from spin-down. The prediction according to the measured v and v in Table [2] is shown as a dotted line. Middle-left 
panel: The effect of the constant slow-down is subtracted from the rotational frequency. Especially the second glitch can clearly be seen 
as a deviation from the parabola-shape. The latter indicates a significant and stable braking-index and the prediction according to the 
measured v is shown as a dotted line. Bottom-left panel: The difference between the measured spin frequency and the contributions from 
V, V and ii. The solid and dotted lines show the prediction according to the glitch model including/excluding glitch recovery. Top-right 
panel: The measured spin-frequency derivative as a function of time. The dotted line indicates the i/ and u contribution. Middle-right 
panel: The difference between the data and the dotted line of the top-right panel. The solid line indicates the prediction according to 
the glitch model including glitch recovery. Bottom-right panel: This plot is identical to the middle-right panel, but slightly zoomed in. It 
can clearly be seen that the exponential glitch-recovery is over relaxing resulting in a positive step in spin- frequency derivative, i.e. the 
pulsar spins down slower than what could be expected from its braking index. 



be done by calculating a TOA for both peaks and by de- 
termining which corresponding residual follows the general 
trend of the residuals closest. However, the double-peaked 
profile was detected in the first observation after the 2007 
glitch, hence during the fast part of the recovery process. 
This makes it impossible to distinguish between the resid- 
uals corresponding to the two peaks because both residuals 
can be explained equally well by slightly adjusting the glitch 
parameters. The timing solution presented in Table [5] there- 
fore makes use of the alignment discussed in Sect. 15.11 

6.1 Glitch properties 

Fig. [To] shows the rotational evolution of PSR J1119-6127 
over more than 12 years. Its spin parameters as a function of 
time were measured using short stretches of data. For each 
TOA we took all TOAs within a range of 75 days and mini- 
mized the residuals using a timing model that only included 
(/)o, V and z> as free parameters and to as a fixed parameter. 
Fig. [To] only includes points which are based on more than 
three TOAs. 



The top left panel of Fig.[TO]shows the frequency against 
time, which decreases steadily according to the measured z>. 
However, the spin-down rate varies slightly with time and 
the effects of /> are only visible when subtracting the efi'ects 
of V first, as shown in the middle- left panel. The almost 
perfect parabola-shape implies a rather constant and sta- 
ble second spin- frequency derivative, from which the brak- 
ing index can be estimated. In this panel it is also possi- 
ble to see the effects of the two large glitches mentioned 
above. Both glitches exhibit recoveries towards the pre-glitch 
trend, as commonly ob served in many young pulsar glitches 
(jShemar fc Lvnelll996l l. However, it can be observed that 
especially after the 2007 glitch the rotation evolves towards 
a new state, involving a lower spin-down rate l^l and there- 
fore showing a larger spin-frequency than the pulsar would 
have if it continued evolving as before the glitch (see the 
bottom- left panel). This particular behaviour is best appre- 
ciated by looking at the evolution of i>, as shown in the 
top-right panel. The spin-down rate evolves almost linearly, 
as expected for a stable i>, only interrupted by the two ma- 
jor glitches and their recoveries. After each glitch there is a 
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Table 2. Rotational parameters for PSR J1119-6127. 



Parameter 


Value 


Epoch (MJD) 


52109.85 


V (Hz) 


2.4512027814(5) 


u (10-15 Hzs-i) 


-24155.07(1) 


u (10-24 Hzs-2) 


638.9(4) 


00 


37.61(2) 


DM (cm-3 pc) 


713 


n 


2.684(2) 


MJD range 


50850-55364 


2004 glitch parameters 


Glitch epoch 


53290 




4.8(3) 


Ai/p (^tHz) 


-0.005(13) 


A;>p (10-15 Hzs-i) 


6.4(3) 


Aj/rf (mHz) 


0.8(1) 


rd (days) 


63(8) 


2007 glitch parameters 


Glitch epoch 


54240 


A(/) 


22(3) 


Afp (mHz) 


-1.0(1) 


Ai>p (10-15 Hzs-i) 


30(1) 


Ai.(i' (/.Hz) 


5.1(1) 


T^^-' (days) 


214(7) 


/\uf^ (/.Hz) 


9(2) 


r^^' (days) 


23(3) 



jump in spin-down rate followed by an exponential recovery 
to pre-glitch values. However, as indicated by the positive 
permanent jump AOp (Table [2]), after the recovery of the 
2007 glitch the spin-down rate is smaller than the projected 
pre-glitch value by almost 30±1 x Hzs-^ (see the last 

two panels at the right hand side of Fig. 1101) . The recovery 
of the 2004 glitch also presents this behaviour, as indicated 
by its positive A!>p = 6.4 ± 0.3 x lO-^^ Hzs-^. 

This post-glitch behaviour is not normal among the rest 
of the pulsar population and it has been clear ly observed 
only for the RRAT PSR J1819-1458 l|Lvne et ah 2009). Nor- 
mally, if a glitch recovery is detected, Ai>p ^ 0, i.e. the spin- 
down rate either goes back to values similar to those before 
th e glitch or stays a t higher rates for a longer time (see Fig. 6 
in iLvne et al.ll2009D . 

The parameter Q = Az^d/Az^g is used to quantify the 
degre e of recovery of the spin frequency (e.g. lShemar fc Lvnd 
Il996t ) and it is observed that p ulsars exhibit Q < 1 for most 
glitches (e.g. iLvne et al.ll200ol ). Because some of the glitch 
parameters depend on the glitch epoch, the uncertainty on 
Q is estimated by obtaining another 4 solutions (in addition 
to the one presented in Table [2]) by varying the glitch epochs 
in the range defined by the TOAs just before and after each 
glitch. The Q parameters obtained are both greater than 
one: Q = 1.01 ± 0.01 for the 2004 event and Q = 1.1 ± 0.1 
for the 2007 one. These anomalous values are due to the 
negative permanent frequency jumps Ai/p observed for both 
glitches. 

A new parameter Q' = ^i>dl ^i>g can be introduced to 
quantify the recovery of the spin-down rate. By making use 
of the relation Ai>d = —Aud/rd, it can be re-written as 

Q'= , . . (8) 

Because Ai/p is usually a negative quantity, normal glitches 



satisfy Q' < 1. The Q' values of the 2004 and 2007 glitches 
are 1.046 ± 0.004 and 1.01 ± 0.02 respectively. As with Q, 
uncertainties are derived by using the 5 different timing so- 
lutions (see above). It is worth noticing that if Avp > 0, as is 
the case for these two glitches, then Q' will be always greater 
than unity. The amount of overpassing reflects the relative 
size of the recovered jump Avd respect to the instantaneous 
jump Aiyg . Consequently, even though the anomalous recov- 
ery is more dramatic for the 2007 glitch (large Avp), Q' is 
larger for the 2004 glitch. 

We note that both Q and Q' are strongly dependant 
upon the exponential fit, including the glitch epoch. Changes 
in the model or the addition of more post-glitch data could 
therefore modify their values considerably. Nonetheless, the 
fact they are greater than unity is significant (despite their 
relatively large uncertainties), since it depends purely upon 
the sign of the permanent jumps Aup or Avp which are 
measured with much higher confidence. In other words, the 
anomalous behaviour of these two glitches is already well 
described by the signs of the pe r mane nt jumps. 

We note that lCamilo et ahl (|2000l ) report a small glitch 
that occurred in 1999 at MJD ~51400. This glitch is much 
smaller {Avg = 0.011 ± 0.001 /tHz) than the glitches we 
report here. Although there is clearly a timing deviation at 
that epoch, Figs. [9land ll0l do not show convincing evidence 
that this glitch is significantly different from other timing 
irregularities. Indeed, there was no need to include this glitch 
in the long term timing solution (Table [2]). 

6.2 Long-term spin evolution 

The long term spin evolution of pulsars can be characterised 
by the so-called braking index n, defined by the power law 
!> = —kf", where fc is a positive constant. With this def- 
inition the braking index can be estimated from the ob- 
servable rotational frequency and its first two derivatives, 
since differentiation of the power law gives n — vvjiP' . 
ferent braking mechanisms w ill produce different values of n 
(Blandfor d fc Romanll988l '). and it is found that most mea- 
sured braking indices are less tha n 3, the canonical value f or 
pure electromagnetic torque (e.g. [Livingstone et al.ll2007l '). 

As it happens with most known very young pulsars, the 
variation of the spin-down rate (i>) is signifi cant and sta- 
ble, a nd therefore relatively easy to measure. ICamilo et al.l 
l|2000t ) measured this variation for PSR J1119-6127 and 
reported a braking index n — 2.91 ± 0.05. Using the tim- 
ing solution presented in Table [2] we obtain a braking index 
n = 2.684 ± 0.002, the error calculated by propagating the 
individual errors on the frequency and its de rivatives. While 
the value reported by ICamilo et al.l |20oO) was measured 
using a data span of about 1.2 yr, this new estimate rep- 
resents the spin behaviour of the pulsar for more than 12 
years, hence the improvement of the uncertainty. The two 
large glitches interrupted the spin-down evolution carried on 
during the first half of the data (see Fig llOp , and it is not 
obvious how they will affect the long term spin evolution. 
If only data until the first large glitch is used, the result is 
n — 2.686 ±0.002, well contained between the error margins 
of the value given above. This similarity makes us confident 
that the glitches are modeled in a physically sensible way. 

We measure a braking index n = 2.684 ± 0.002, smaller 
than the predicted value for pure dipolar electromagnetic 
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braking, which imphes the presence of extra torques. The 
long-term spin evolution of this pulsar is however unclear, 
and may not be solely dictated by the measured braking 
index. If glitches like the one reported here are normal, the 
effective long-term evolution would have to be described by 
a larger braking index. In that case, the quoted value would 
better represent the inter-glitch, rather than the long-term 
spin evolution. 



7 DISCUSSION 

7.1 The glitch-induced identity changes of PSR 
J1119-6127 

One of the most intriguing properties of PSR J1119-6127 
is that the radio emission of its distinct profile components 
can be classified differently. Its profile is most of the time 
single peaked and in no way remarkable. Only once in more 
than 12 years of data a very different pulse profile was ob- 
served with a second transient peak appearing next to the 
"normal" leading peak. The switching timescale is therefore 
more similar to those of the intermittent pulsars than the 
shorter timescale mode changes observed in many other pul- 
sars. The leading component in the transient double-peaked 
mode is usually active at radio wavelengths and it is the 
trailing peak which has emerged. This identification is made 
by comparing the PA curve in the normal and the double- 
peaked mode. In addition, the circular polarization and the 
profile shape show a similar correlation. 

Radio emission which can be put into a third category 
was discovered by analysing the individual pulses of this 
source. In a single observation made at a wavelength of 20 
cm four strong pulses leading both the "normal" and "in- 
termittent" component were detected. These pulses are spo- 
radic enough to make the signal disappear in the noise when 
adding the emission of all rotations during the observation 
together, thereby making this a "RRAT-like" component. 
At a wavelength of 10 cm a second component showing ex- 
tremely erratic emission was detected trailing the other pro- 
file components, thereby making these erratic components 
flank the other emission. 

The "identity changes" of PSR J1119-6127 are clearly 
glitch-induced. Not only it has been demonstrated that the 
rate of occurrence of the RRAT-like pulses at a wavelength of 
20 cm must have been significantly higher around the epoch 
of their detection, but more importantly, all abnormal be- 
haviour is observed directly after the large amplitude glitch. 
The transient component was detected in 2007 June 17 (once 
in 49.8 hours of data), which is the first observation after the 
glitch. The leading RRAT-like component was detected at a 
wavelength of 20 cm in 2007 August 20 (only four pulses out 
of the 246,892 inspected stellar rotations in total). This ob- 
servation was the first long individual pulse recording after 
the glitch at a wavelength of 20 cm (3211 pulses long). Three 
other 20 cm observations were made between this observa- 
tion and the glitch epoch. However, with a total of only 945 
recorded stellar rotations and the relatively low event rate 
it is not surprizing the RRAT-like pulses were only detected 
during the longer observation. Finally, the trailing and lead- 
ing erratic components were detected in 2007 July 23 at a 
wavelength of 10 cm, which is the first 10 cm observation 



after the glitch. This makes it likely that erratic RRAT-like 
pulses were emitted throughout the time between the glitch 
epoch and their detection at 20 cm. 

Because of the separation between the observing ses- 
sions, the glitch epoch cannot be accurately determined. As 
far as we can tell, the abnormal emission properties of PSR 
J1119-6127 started directly after the glitch and continued 
for 78±24 days. This corresponds to 4±1 times the timescale 
of the fast part of the recovery or 0.4±0.1 times the timescale 
of the slow part of the recovery. This is the first time that a 
glitch, or the post-glitch recovery process, is observed to in- 
fluence the radio emission process of a normal (non-RRAT) 
pulsar. 

It is currently not possible to know for sure how dif- 
ferent the outer erratic components are compared to the 
intermittent second peak. Both phenomenon are observed 
around the same epoch and appear to be transient events 
and could therefore represent different phases of the same 
phenomenon. The additional components probably repre- 
sent regions on the polar cap which usually fail in producing 
the positron-electron pairs required for the emission mech- 
anism to work. One could speculate that a large-amplitude 
glitch for some reason changes the magnetosphere or stellar 
surface such that these regions switch-on, either completely 
or in a flickering manner. Possibly, for all large glitches the 
steady emission is followed by a flickering state, before the 
physical conditions relax back to the pre-glitch state and the 
components switch off again. It must be noted that the over- 
all profile symmetry with the more steady components in the 
centre and the more flickering components at the edges of 
the beam suggests that the different components do have 
distinct emission properties rather than we observing the 
different components in different phases of the same cycle. 

It is not clear why the occurrence of a glitch and the 
changes in the radio emission are related to each other. If 
the glitch is triggered solely by the conditions inside the star, 
than how is the energy involved in the glitch transported to 
the magnetosphere to induce the changes in the emission 
properties? This problem could suggest that there is an ex- 
ternal event (a change in the magnetosphere) that triggers 
both the glitch and the changes in the emission properties. 
In the next subsection we will argue that the glitches might 
be triggered by magnetic stresses on the crust. 

7.2 The glitches of PSR J1119-6127 

7.2.1 Companson with other pulsars 

Detected glitch frequency jumps Avg are in the range 0.001- 
50 /iHz among the pulsar population, im plying fractional 
chan g es /S.v„/v between 10^^^ and 10~^ IjShemar fc Lvng 
1 19961 : IWang et al]|200Gl ). Young pulsars Uke the Crab and 
PSR J1119— 6127 normally exhibit smaller glitches, with 
l^v/v < 300 X 10"^. However, the 2007 glitch, with Ai/g/i/ ~ 
5400 X 10~®, is as big as the largest glitches detected in the 
rotation of the Vela pulsar and many other young pulsars 
with characteristic ages between 10 and 100 kyr. 

Exponential glitch recoveries have been observed be- 
fore for other pulsars and some of them, like PSR J1119- 
6127, have required more than one time constant. For in- 
stance, more than one timescale was n ecessary to mode l 
the very large glitch in PSR B2334-I-61 jYuan et aLlbOld ) 
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Table 3. Recovery parameters for PSR J1119— 6127 and other pulsars. Uncertainties for the last quoted digit are given in parenthesis. 
Note that the fact that Q and Q' of PSR J1119— 6127 are larger than one is significant, despite the large error bars on the values 
themselves (see text). 



PSR 


MJD 








Q 


Q' 


Refs. 




[days] 


(10-6) 




[10-1^ Hzs-i] 








J1119-6127 


54240 


5.4(8) 


-1.0(1) 


30(1) 


1.1(1) 


1.01(2) 




J1819-1458 


53924.8 


0.699(4) 


0.138(1) 


0.79(1) 


0.159(4) 


1.7(5) 


1 


J1846-0258 


53883 


4(1) 


-95(1) 


-274(13) 


9(3) 


0.973(2) 


2, 3 


B0355+54 


46497 


4.368(2) 


27.8773(1) 


-0.5(1) 


0.00117(4) 


0.14(3) 


4 


B1800-21 


48245 


4.07(2) 


30.07(1) 


-40(1) 


0.0137(3) 


0.46(2) 


4 


B2334+61 


53615 


20.579(1) 


41.35(3) 


-8.68(2) 


0.0075(1) 


0.92(1) 


5 



Refs.: l.lLvne et al.l ||2009|'I: 2. [Livingstone et all a2010fl: 3. [Livinestone et al.n2006i1: 4.iShemar fc Lvn6l"1l99a): 5.lYuan et al.l ||20100 



and glitches in the Vela pulsar (|McCulloch et al] Il990l : 
lAlpar et al.lll993l ). 

A comparison with PSR J1846-0258 is particularly in- 
teresting because its spin parameters (jc — 0.8 kyr and a 
magnetic field of ~ 4 — 5 x 10^'' G) are very similar to that 
of PSR J1119-6127. Until now, this rotationally powered ra- 
dio quiet pulsar was the only very young pulsars for which a 
glitch of this magnitude had been detected. As observed in 
PSR J1119— 6127, the large glitch on that pulsar presented 
a very efficient recovery, with Q — 9 ± 3, thanks to a large 
negative perm anent frequency jump Aup — —95 ± 1 /xHz 
l|Livingstone e t al. 2010). In addition to this anomalous re- 
covery, the pulsar exhibited a radiative outburst episode in 
X-rays resembling magnetar activity, and coincident with 
the l arge glitch llGavriil et all l2008i : iKumar fc Safi-Harbl 
l2008l : lKuiper fc Hermsenll2009l ). 

Glitches are believed to represent changes in the con- 
figuration of the star itself and in general are not associated 
to external factors. No changes on pulse shape or radiative 
behaviour have been reported for normal rotation powered 
pulsars, with the only exception being PSR J1846— 0258. 
One glitch in the RRAT PSR J1819-1458 also exhibited 
an anomalous spin-down recove ry, with Ai>p = .79 ± 1 x 
10~^^ Hzs"^ and Q' = 1.7±0.5. [Lvne et al.l (|2009l ) reported 
a possible increase in the pulsed emission rate accompanied 
by an increase of radio fluxes of the pulses associated with 
this glitch. 

Despite the fact that these three neutron stars have 
shown anomalous glitch recoveries, they have somewhat dif- 
ferent properties. Table [3] displays the values Aup, A/>p, Q 
and Q' for the largest glitches in each of these three objects. 
Additionally, to compare with some other normal glitches, 
we also show the same para meters for the ver y large glitch in 
the young PSR B2334-I-61 jYuan et al.ll20l d). a large glitch 
in the Vela-like pulsar PSR B1800— 21 and a large glitch 
in the relatively old PSR B0355-I-54. Values were estimated 
from published fitted parameters and references are given 
in the last column. While the glitch in PSR J1119-6127 
has both Q and Q' greater than unity, the RRAT exhibits 
Q' > 1, but Q < 1. On the other hand, the glitch in PSR 
J1846— 0258 presents the opposite situation, with Q > 1 but 
Q' < 1. The other two normal young pulsars present nor- 
mal Q < 1 and Q' < 1, owing to normal permanent jumps 
Aup > and Az>p < 0. 

The emission changes observed in PSR J1119— 6127 fol- 
lowing the 2007 glitch (Sect. [3] and [4]) and the anomalous 
recovery the glitch presented are suggestive of a single con- 



nected event, as observed in PSR J1846-0258 and in RRAT 
PSR J1819— 1458. Both these last objects are high magnetic 
field neutrons stars, just as PSR J1119— 6127. 



7.2.2 Glitch models 

The most simple and intuitive model of a glitch is a 
starq uake produced by sudden cha nges of the star shape 
(Bav m et al. I ll969l : lRuderman|[l969l ). These changes can be 
produced regularly as the star slows down and its shape 
evolves from an oblate configuration towards a more spher- 
ical one. However, the high glitch activity of some pulsars 
(like V ela or PSR B1737— 30) ca nnot be accounted for by this 
model (|McKenna fc Lvnelll99d ). A second model solved this 
problem by involving the physics of the interior of the star. 
In this model glitches are produced by a sudden transference 
of angular momentum from the inner neutron superfluid to 
the solid star crust. Any rotating superfluid slows down via 
migration away of the rotation axis of quantised vortices; 
carriers of the superfluid angular momentum. If the migra- 
tion is stopped the superfluid retains its angular momentum 
and does not slow down. This model assumes that vortices 
often pin to different places of the inner crust due to in- 
teraction with the neutron superflui d associated with the 
crustal ion lattice l|Alpar et al.l Il984l ) . In this way the in- 
ner superfluid keeps rotating at at higher rate and behaves 
as an angular momentum reservoir. The collective unpin- 
ning of many vortices may occur once stresses have built up 
due to the increasing differential rotation between the in- 
ner superfluid and the crust, produ cing the sudden spin up 
of the crust we observe as a glitch (jAnderson fc Itohlll975l : 
iMelatos fc Warszawskil 12009). 

Another model, proposed by iRuderman et al.l l|l998h . 
states that vortex migration would push magnetic flux tubes 
away, producing stresses in the rigid crust that may be re- 
leased through quake- like events. In addition to produce pos- 
sible changes of the moment of inertia of the star, these re- 
configurations could trigger the unpinning of a large number 
of vortices, producing a glitch. 

Glitches in magnetars present a broad range of prop- 
erties. Among them, the coincidence between radiative out- 
burst and some glitches is the most relevant, but also the un- 
usual recove ries and sometim es large post-glitch spin-down 
increments (|Dib et al.l |2008| ). These characteristics, which 
are not seen in normal radio pulsars, have induced the 
hypothesis that crust rearrangements caused by magnetic 
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field decay or reconfigurations may tri gger these gUtches 
l|Dib et al.ll2007l : iLivingstone et al.ll2010l ) . 



7.2.3 Interpretation of glitch recoveries in terms of the 
pinning -unpinning model of glitches 

The pinning-unpinning model involves the existence of two 
main components: a resistive one, composed of all vortices 
which are continuously unpinning and pinning, and a ca- 
pacitive one, containing those vortices that unpin only dur- 
ing a glitch (|Alpar et al.lll996l ). The resistive component is 
continuously transferring angular momentum to the crust, 
as vortices unpin due to thermal fluctuations or quantum 
tunnelling (|Alpar et al.lll984l : iMelatos fc Warszawskill2009^ . 
The capacitive component does not participate in the regu- 
lar spin down of the star until a glitch occurs, and the star 
may contain several capacitive zones. 

During a glitch the crust spins up and the rotational 
lag with the superfluid decreases, causing the effects of the 
resistive component to decrease l|Wong et alj bOOll ) . As a 
consequence, the spin-down rate increases due to the action 
of the always present external torque. This explains the ini- 
tial spin-down jump at the glitch epoch and the following 
relaxation towards the equilibrium configuration, as vortices 
from the resistive component continue to pin and unpin. 

Once the recovery is finished one would expect the star 
to follow the same spin-down rate as before the glitch. How- 
ever, most young pulsars for which a recovery is observed 
present positive permanent jumps in frequency derivative. 
For the Crab pulsar, it has been proposed that this change 
in spin-down is caused by the creation of vortex trap zones, 
where many vortices pin and keep out of the regular slow 
down (jAlpar et al.lll996l : IWong et al.ll200ll ). The pinning of 
a large number of vortices decreases the effective moment 
of inertia over which the external torque is acting, hence in- 
creasing the observed spin-down rate. In other words, during 
the glitch many vortices unpin, producing the observed spin 
up, but also a larger number of vortices pin in trap zones. In 
this context, Q values less but similar to one are interpreted 
as a capacitive behaviour, where the pulsar is building an 
angula j moment reservoir l|Wong et al.l l20oj] '). lAlpar et al.l 
l|l996f ) state that young pulsars like the Crab are still cre- 
ating their capacitive components, which later will produce 
the larger glitch activity observed in slightly older pulsars 
like Vela. Additionally, they propose that the smaller and 
less frequent glitches in the Crab pulsar are likely to be trig- 
gered by starquakes, produced by the rapid spin-down of the 
star. 

Following these ideas, and because of their similar age, 
the glitch activity in PSR J1119-6127 should present sim- 
ilar characteristics as the Crab. However, we have shown 
that the two last glitches of this pulsar present anomalous 
characteristics, specially on their spin-down rate recoveries. 
If the number of vortices that unpinned was considerably 
larger than the number of vortices that pinned (during the 
creation of capacitive zones), then the effective moment of 
inertia after the glitch would be larger than before, and 
the spin-down rate smaller. This would produce a post- 
glitch spin-down evolution like the one observed in Fig. 1101 
What could cause this different behaviour is not clear, but 
the high magnetic field strength is something that PSRs 
J1119-6127, J1846-0258 and RRAT J1819-1458 have in 



common. Maybe, these glitches have been triggered by mag- 
netic stresses, which could deform the crust, thereby chang- 
ing the moment of inertia, or release pinned vortices. Either 
way, such a glitch would not be part of the regular glitch 
activity, and those processes following a normal glitch may 
not necessarily occur. An unusually large number of vortices 
may have unpinned in PSR J1119— 6127 (considering the low 
glitch activity of pulsars of its age), producing a very large 
spin up. But, because the origin of the glitch was different, 
the unpinning of vortices overcame any possible vortex pin- 
ning, making up a larger final effective moment of inertia 
and therefore a smaller spin-down rate. 



7.3 Comparison with other neutron star 
populations 

7. 3. 1 Comparison with other young pulsars 

The pulse profile of PSR J1119-6127 in its double-peaked 
mode (see middle panel of Fig. [S} has striking similari- 
ties with a group of young pulsars with a very distinct 
pulse shape. These energetic wide-beam pulsars, such as 
PSRs J1015-5719, J1105-6107 and J1420-6048, are char- 
acterized by their relatively wide double-peaked profile with 
steep inner edges and a hig h degree of linear polarization 
([Johnston fc Weisberell2006l ). The similarity of this type of 
profiles with that of PSR J1119-6127 therefore suggests that 
the fiducial plane corresponds to the centre of the double- 
peaked profile. Having the magnetic axis in the middle of 
the double-peaked profile makes it more natural that the 
two peaks have the same amplitude when they are visible. 
This scenario (case IV of the discussed geometries in Table 
[1]) was also preferred after comparing the pulse profiles at 
the different observing frequencies. 

The main difference between this pulsar and other pul- 
sars of the same class of profiles might be that the PA curve 
is not as flat under the flrst of the two peaks in the double- 
peaked mode. The estimated emission height of ~500 km 
is som ewhat larger than that of t he other pulsars in this 
group (|johnston fc Weisberd [20061 derived emission heights 
of 380, 110 and 175 km for PSRs J1015-5719, J1105-6107 
and J1420-6048 respectively). The relatively wide observed 
pulse proflle is attributed to a relatively aligned geometry 
(a ~ 17° — 30°), or possibly smaller if only part of the open 
field line region is active. 

There are a few young energetic pulsars which exhibit 
relatively strong individual pulses, although they are not 
glitch related. The strong pulses at the leading edge o f the 
pulse profile of the Vela pulsar ('Johnston et alj 120011) and 
at th e trailing edge of PSR B1706-44 ( Johnston fc Romanil 
l2002t ) were dubbed "giant micro-pulses" . These giant micro- 
pulses have large peak-fiux densities, but their mean flux 
densities are not as extreme as t he giant pulses observ ed for 
instance for the Crab pulsar (e.g. lLundgren et al.lll995l ). Un- 
like normal emission, the pulse-energy distribution of both 
giant pulses and giant micro-pulses are best described by a 
power law and are therefore thought to have a somewhat dif- 
ferent origin. Unfortunately, not enough strong pulses were 
detected for PSR J1119-6127 to quantify their pulse en- 
ergy distribution. Nevertheless, there are characteristics the 
strong pulses of PSR J1119-6127 do share with the giant 
micro-pulses: they all occur at the edges of the pulse proflle 
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and are broader than giant pulses (e.g. lHankins et al]|2003h . 
Another characteristic that Vela and PSR B1706-44 have in 
com mon is that they are both detected at 7-ray energies 
fe.g. lAbdo et al.ll201(]| ). The E of PSR J1119-6127 suggests 
that this pulsar might be a 7-ray source as well. 



7.3.2 Comparison with the RRATs 

It can be argued that if the line of sight to this pulsar were 
different, one would only be able to detect this pulsar via its 
individual pulses, hence it would have been classified as a 
RRAT. This would suggest that some of the RRATs are very 
similar in the sense that they may emit normal radio pulses, 
but not in the direction of Earth. There are more analogies 
one can draw between PSR Jf lf9-6f 27 and the RRAT pop- 
ulation. The unusual post-glitch over-relaxation observed for 
PSR J1119-6127 is not normal among the rest of the pulsar 
populati on, but it is simil ar to that observed for one of the 
RRATs (|Lvne et al.ll2009l '). Also, three of the seven RRATs 
with a measured P have a Bs > 10^'^ G (|McLaughlin et al] 
I2OO9I ) and have therefore field strengths comparable with 
that of PSR J1119-6127. In addition it must be noted that 
the evolutionary track of PSR J1119-6127 in the P - P dia- 
gram may lead it towards the RRAT population rather than 
the bulk of the normal pulsars as is suggested by its braking 
index. 

Having said this, a difference between PSR J1119-6127 
and the RRATs is that in general there is no evidence that 
the event rate of RRAT pulses is time-dependent, certainly 
not at the level as observed for PSR J1119-6127. One of the 
RRATs is argued to hav e an increase in th e RRAT-activity 
associated with a glitch (|Lvne et al.ll2009l ). although unlike 
PSR J1119-6127 it is always observed to be RRAT-like. If 
more sources like PSR J1119-6127 exist, there might be a 
population of neutron stars which emit pulses extremely spo- 
radically after a glitch, possibly only some clustered pulses 
every so many years. This kind of s ources, like "Lorimer 
burst" type of events (|Lorimer et al.1 ,2007). would be hard 
to detect. It will probably require next-generation telescopes 
(such as the LOFAR, ASKAP, MeerKAT or the SKA), which 
can afford long dwell times because of their large field-of- 
view, to discover this type of object. 

PSR J1119-6127 shows erratic emission on both sides 
of the profile cone. This might be related to the general 
observation of the longitude resolved modulation index be- 
ing higher at the edges of profile components (see e.g. 
iTavlor et al1ll975l : IWeltevrede et aLll2006l ). This "edge ef- 
fect" suggests that the locations of the profile peaks cor- 
respond to field lines at which the emission mechanism is 
stable. The edge effect is clearly obse rved for the Vela pul- 
sar, PSR B1046-58 and B1706-44 djohnston et al.1 12001I: 
Johnston fc Romanil [20o3) . PSR B1133-I-16 (|Kramer et all 



20031 ). and PSR B0656-I-14. The latter pulsar shows the 



strongest pulses in the centre of the profile (up to 116 times 
the average pulse intensity), but the emission on the leading 
edge can be up to a factor 2000 b righter than the averag e 
emission at those pulse longitudes (| Weltevrede et aLll2006l ). 
Some of the RRATs might therefore be pulsars for which the 
line of sight only intersects this unstable part of the emission 



7.3.3 Intermittent pulsars 

The double-peaked profile is only observed once in 12 years 
of data. This ti mescale is more si milar to those of the inter- 
mittent pulsars (|Lvne et al.ll2O10l ) than the shorter timescale 
mode changes obse rved in many other pulsars (e.g. ,Backerl 
I1970I : [Rankin|[l986l ). Apart from the similarity in timescale, 
the most obvious connection between PSR J1119-6127 and 
the intermittent pulsars is that changes in profile shape are 
linked to the timing properties of the star. However, the 
physical mechanism of this link are likely to be different. 
In the case of the intermittent pulsars it is believed that 
changes in the magnetosphere affect the pulsar rotation (the 
value of v) , while for PSR J1119-6127 we argue that it is pos- 
sibly an external change in the magnetosphere that triggers 
the glitch. In any case, both the intermittent pulsars and 
PSR J1119-6127 show that there is an intimate connection 
between the neutron star interior and the magnetosphere. 



7.3.4 Comparison with PSR B0656+U 

PSR J1119-6127 shares an important feature with PSR 
B0656-I-14. Both pulsars were discovered using standard 
periodicity search techniques, but subsequently found to 
exhibit RRAT-like emission. Both pulsars therefore pro- 
vide a link between the "normal" pulsars and the RRAT- 
like pulsars. Pul sed X-rays have bee n observed for both 
PSR J1119-6127 (iGonzalez et al.|[2005l') and PSR B0656-H4 
(|C6rdova et al.l Il989l : iDe Luca et al.l I2OO5I ). In the radio 
band they both have a more or less triangular-shaped pulse 
profile which covers only a part of the open field line region 
(in the case of PSR B0656-I-14 there is a trailing shoulder, 
while PSR J1119-6127 has the intermittent trailing peak). 
They both have a high degree of linear polarization at 1.4 
GHz and it will be interesting to see if PSR J1119-6127 be- 
comes depolarized at higher frequencies, like PSR B0656+14 
(Johnston ct al. 20QS). 

There are also some important differences between 
PSRs J1119-6127 and B0656-hl4. PSR B0656-hl4 is a 
middle-aged pulsar located near the bulk of pulsars in the P- 
P diagram, while PSR J1119-6127 is very young and has an 
extremely strong magnetic field. The RRAT-like emission of 
PSR J1119-6127 is found in an separate profile component 
(more like giant micro- pulses) , while the strongest pulses of 
PSR B0656+14 are found at the same pulse longitude as the 
normal emission. Therefore the RRAT-like emission of PSR 
B0656-I-14 is interpreted as being the tail of an extended 
pulse-energy distribution ([Weltevrede et al]|2006l ). while in 
the case of PSR J1119-6127 and the giant micro-pulses the 
RRAT-like emission is more distinct from the normal emis- 
sion. Also, the emission of PSR B0656-I-14 was found to 
be more extreme at low observing frequencies, while for 
PSR J1119-6127 strong pulses possibly occur more often at 
higher observing frequencies. Therefore the physics of the 
strong pulses observed for PSR B0656-I-14 might be more 
related to th ose observed for older pu lsars at very low fre- 
quencies (e.g. lErshov fc Kuzminll2005D . In addition to these 
differences the strong pulses of PSR B0656-f 14 are not re- 
lated to glitches. 
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7.5.5 Comparison with magnetars 

Because of the high magnetic field strength of PSR J1119- 
6127 there is an obvious direct link with the magnetars. 
Links between magnetars and rotation-powered pulsars have 
been suggested before as for instance PSR J1846-0258, a 
young (radio-quiet) pulsar with a similar magnetic field 
strength to that of PSR J 1119-6127, h as shown magnetar- 
like X -ray bursts (|Gavriil et all l2008l : iKumar fc Safi-Harbl 
20081). In addition, an over-recovery (Q = 8.7 ± 2.5; 



Livingstone et al]|2010l ) has been claimed for that source, 
as well as for magn etar AXP 4U 0142+61 {Q = 1.07 ±0.02; 
iGavriil et ai]|2009D . 

T hree of the magne t ars ar e detected as radio pul- 
sars (jCamilo et all |2006| . l2007l : iLevin et all l2010h and 
they show characteristics similar to PSRs J1119-6127 and 
B0656-I-14. In part icular all three have a high degree of 
linear polarization jCamilo et al.l [20071 : iKramer et al.ll2007l : 



ICamilo et aIll200S ) and individual pulses from PSR J1810- 
197 l|Servlak et al.P 2009') and J 1623-4950 (|Levin et al.ll2010l ') 



show spiky emission sim ilar to that seen in PSR B0656-I-14 
l|Weltevrede et al.ll2006l ). Curiously, evidence seems to sug- 
gest th at the magnetars arc close to being aligned ro- 
tators (|Camilo et all |2007| . 2008) and this might also be 
the case for PSR 1119-61 27 and possibly PSR B06564-14 
l|Everett fc Weisberd [200 ll ). Could it be that the aligned 
geometry is responsible for some of these attributes? We 
note als o that the 7-ra y profile of PSR B0656-I-14 is pecu- 
liar (W eltevrede et al.l l2010!') and we await the detection of 
PSR J1119-6127 with the Fermi satellite with interest. 



7.4 What are RRATs? 

The inconsistency between the observed supernova rate and 
the size of the RRAT population sugges ts there are links 
between the different neutron star classes ijKeane fc Kramer! 
|2008| ). Indeed, it has been argued that one of the RRATs, 
PSR J1819-1458, might be evolved from the magnetars 
based on the unusual P change after its glitch (jLvne et al.l 
I2OO9I ). Because of the movement of the Crab pulsar in 
the P — P diagram the youn g pulsars an d the magnetars 
appear to be linked as well (jLvnel [20041 ) . The view that 
the RRATs are an evolutionary stage rather than a sepa- 
rate class of neutron stars is further strengthened by the 
discovery of "missing links": pulsars which can be classi- 
fied both as pulsars or as RRATs. Some RRATs will turn 
out to be pulsars with exten ded pulse energy distrib utions 
(like that of PSR B0656-H4; IWeltevrede et al.ll2006l ). oth- 
ers might have a "normal" profile component (like PSR 
J1119-6127) which happen to be missed by our line of sight, 
others might be pulsars with extreme null lengths (see e.g. 
iBurke-Spolaor fc Bailesll20ld ) and some might be magnetars 
with more extreme spik y radio emission than XTE J1810- 
197 (jServlak et al.ll2009l ). 

With the above in mind, a valid question to ask one- 
self is what you are left with when you remove the pulsars 
that fall in the above categories from the RRAT population. 
Will there be any "true" RRATs left? Or are RRATs sim- 
ply the combination of the extreme ends of the other classes 
of neutron stars? Because RRATs do not form a separate 
island in the P — P diagram and because it proves to be 
difficult to define the RRATs based on their physical prop- 



erties (see e.g. iKeane et al.ll2O10l ). it is tempting to apply 
Occam's razor and assume that the RRATs do not form a 
distinct population of pulsars. 

A problem with the term RRAT is that it is not clearly 
defined. The definition that RRATs are pulsars which are 
not detectable via periodicity searches is insufficient because 
(as illustrated by PSR B0656-I-14) it depends on the dis- 
tance of the source (and the integration time of the observa- 
tion and the sensitivity of the telescope). By accepting that 
RRATs are not a separate class of pulsars the question of 
how to define a RRAT is irrelevant. More useful would be 
to define how erratic (hence "RRAT-like") the emission of 
a pulsar is. Such a definition would allow a better investi- 
gation of how this type of emission is correlated with other 
neutron star parameters. 

The degree of RRAT-like behaviour can be quanti- 
fied by intermittency ratio (McLaughlin fc Corded l2003l : 



iDeneva et al]|2009l : [k eane et al.ll2010t ). which compares the 
signal-to-noise ratio of the detection of the pulsar in a 
single-pulse s earch with that of a p eriodicity search, or the 
R-parameter (jjohnston et al.ll200ll ). which is a normalized 
measure of the brightest observed pulse. Both measures have 
the problem that it depends on both the sensitivity of the 
telescope and the duration of the observation. A straight- 
forward way to define the degree of RRAT-like behaviour 
would be to use the modulation index (Eq.[T|), which quan- 
tifies the broadness of an amplitude distribution and is 
observing-system independent. Note that m is usually cal- 
culated as a function of pulse longitude. This makes a big 
difference for especially the magnetar emission, which shows 
many bursts per stellar rotation. 

For RRATs {I) is below the detection threshold and 
therefore the modulation index cannot be calculated (the 
same problem applies for the intermittency ratio or the R- 
parameter). Nevertheless one can determine a lower limit 
on the modulation index by taking zero intensity for the 
undetected pulses and replacing (J) by its maximum allowed 
value, which is the detection threshold for the average pulse 
intensity {I)^^^^^- The resulting limit is 



m ^ 



/ , iVdct N 



1/2 



(9) 



where A'dct is the number of detected pulses. 

The modulation index for selected sources can be found 
in Table |4] One can see that the modulation index of the 
normal emission of PSR J1119-6127 is typical for that of 
radio pulsars. Unfortunately, the lower limit of the modu- 
lation index of the RRAT-like emission of PSR J1119-6127 
is not much higher than that of the normal emission, and 
therefore it is consistent with all types of emission. To im- 
prove this limit one needs to use an instrument which is 
much more sensitive than Parkes, such as the SKA. On the 
other hand, one can see that the modulation index of RRAT 
J1819-1458 must be much higher than that of normal pul- 
sar emission and is possibly more similar to that found for 
magnetar emission or that of giant pulses. Measurements of 
m for more RRATs and other sources with extreme radio 
emission has the potential to differentiate between different 
types of radio emission. 
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Table 4. Table with some typical values of the modulation index 
m measured for different types of emission. 



Source 


Type of emission 


m 


Ref. 


Typical pulsar 


Normal 


0.5 


a 
T 


Crab pulsar 


Giant pulses 


5-8 


PSR B0656+14 


RRAT-like 


~ 5 




PSR J1119-6127 


Normal 


- 0.5 


PSR J1119-6127 


RRAT-like 


^ 0.6 


7 


XTE J1810-197 


Magnetar emission 


1 - 10 


PSR J1819-1458 


RRAT 


^ 6.6 





" Typical minimum value of the lo ngitude resolved m odulation 
index at a wavelength around 21 cm dWeltevrede et al.lliooei 'l and 
m is typically observed to flare up at the profile edges. 
^ Longitude resolved modulation index at the position of the gi- 
ant pulses at 1400 MHz ( Weltevrede ct al. 200j). The actual value 
might be higher if the giant pulse s were not entirely resolved. 

Measured by Weltev rede et all ||2006| ') at a frequency of 327 
MHz in the centre of the profile (where the brightest burst oc- 
curred) . 

This paper. The modulation index flares up at the profile edges. 
^ Modulation index of the RRAT-like component based on the 
single 20-cm observation with the four individual pulses. 

The longitude resolved modulation index is variable in time, 
but in the range of m = 1 — 4 at 1.4 GHz or m = 1 — 10 at 
8.35 GHz jServlak et al.l l2009l. 

3 This is the RRAT with the best determined pulse energy dis- 
tribution. The estimate d limit is based on this d istribution (Sup- 
plementary Figure 2 of iMcLaughlin et al.l [20061 ) in combination 
with the quoted event rate of 229 detected pulses in 13 hours of 
data and the statement that average peak flux density must be 
less than 0.05% that of the strongest detected burst. 



8 CONCLUSIONS 

We found that the young high-magnetic field pulsar J1119- 
6127 exhibits three types of pulsar behaviour. It has an "in- 
termittent" profile peak trailing the peak associated with 
"normal" emission and these profile components are fianked 
by two much more erratic components, which can be argued 
to be "RRAT-like" . Both the erratic components and the in- 
termittent component are observed around the same epoch 
and appear to be transient events. It can therefore not be 
ruled out that they represent different phases of the on- and 
off-switching process, possibly related to regions on the po- 
lar cap which usually fail in producing the positron-electron 
pairs required for the emission mechanism to work. Never- 
theless the overall profile symmetry suggests that the more 
steady components are located in the centre of the beam 
and the more fiickering components at the edges. This is 
consistent with the "edge effect", which states that in gen- 
eral the emission is more erratic at the edges of the profile 
components. Some of the RRATs might therefore be pulsars 
for which the line-of-sight misses the part of the beam with 
more steady emission. 

Both the intermittent and RRAT-like events are very 
rare and are preceded by a large amplitude glitch that 
occurred directly before these phenomenon were observed. 
This is the first time that a glitch, or the post-glitch recovery 
process, is observed to influence the radio emission process 
of a normal (non-RRAT) pulsar. The glitch was followed by 
a post-glitch behaviour which is very unusual for the pulsar 



population as a whole, but it is similar to that observed for 
one of the RRATs (PSR J1819-1458). Aft^er the glitch, the 
spin-down rate relaxed to a smaller value than the projected 
pre-glitch value. We interpret this as an increase of the ef- 
fective moment of inertia produced by an excessive number 
of unpinning vortices. By considering this anomaly and the 
emission changes observed, we believe that the glitch may 
have been triggered by magnetic reconfigurations. 

Although we measure a braking index n = 2.684±0.002, 
glitches like the one reported here make the effective brak- 
ing index larger. Nevertheless, the evolutionary track of 
PSR J1119-6127 in the P - P diagram may lead it to- 
wards the RRAT population rather than the bulk of the 
normal pulsars. The occurrence-rate of RRAT-like pulses 
must have been higher around the glitch epoch (as is claimed 
for J1819-1458), suggesting a link with the glitch-recovery 
process. More precise, the RRAT-like pulses were observed 
during the initial fast recovery (on a timescale of ~ 21 days) 
rather than the later slow recovery (on a timescale of ~ 208 
days). This allows the possibility that RRAT-like activity is 
in general low for young, high-magnetic field pulsars (except 
after a glitch) and increases when pulsars age. 

Like PSR B0656+14, which would be classified as a 
RRAT were it more distant, PSR J1119-6127 provides a link 
between the known neutron star population and the RRATs. 
Additional links could be other young pulsars which show 
"giant micro-pulses", the magnetar population which show 
PSR B0656-|-14-like spiky emission and pulsars with extreme 
long null lengths. It appears that RRATs represents a cor- 
ner in the P — P diagram in which the pulsars have the 
most erratic emission. We argue that the modulation index 
is a good measure of how RRAT-like the emission is. The 
modulation quantifies the broadness of the pulse energy dis- 
tribution and is (unlike for instance the "intermittency ra- 
tio" or the "R-parameter" ) independent of the sensitivity of 
the telescope or observation duration. For RRATs (which 
are not detected in a periodicity search) one can only de- 
rive a lower limit. This limit is estimated to be m ^ 6.6 for 
the RRAT PSR J1819-1458, which although high for pul- 
sar standards is not necessarily higher than that observed 
for magnetars. Additional observations could place stronger 
limits on how high the modulation index of RRAT emission 
is and can potentially be used to distinguish between the 
different sub-classes of the RRATs. 

Not only are the strong individual pulses of PSR J1119- 
6127 reminiscent of those of other young pulsars, but so 
is the shape of its pulse profile in its double-peaked mode. 
The similarity with other "energetic wide-beam pulsar" sug- 
gests that the magnetic axis is located in between the two 
components. The derived emission height for that scenario 
is ~500 km, somewhat higher than that derived for other 
young pulsars, thereby suggesting a relatively aligned ge- 
ometry (q ~ 17° — 30°), or slightly more aligned if only part 
of the open field line region is ever active. 
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